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Abstract 

Consider the polynomial tr {A + tB)™" in t for positive hermitian matrices A and 
B with TO e N. The Bessis-Moussa-Villani conjecture (in the equivalent form of Lieb 
and Seiringer) states that this polynomial has nonnegative coefficients only. We prove 
that they are at least asymptotically positive, for the nontrivial case of AB ^ 0. More 
precisely, we show — once complex-analytically, once combinatorially — that the k-ih 
coefficient is positive for all integer m > toq, where toq depends on A, B and k. 
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1 Introduction 



Some 30 years ago, Bessis, Moussa and Villani (BMV) conjectured [Tj^l that for any hermitian 
n X n matrices A and B, the function 

H{t) := tr exp{A - tB) 

with t G M is the Laplace transform of a positive measure on [0, oo), provided B is positiv^. 
Lieb and Seiringer |12j proved that this statement is equivalent to the assertion that, for 
positive integers m and positive hermitian A and i?, the polynomial 

has nonnegative coefficients only. Here, the Hurwitz product S^.k [S] equals the sum of all 
words in A and i?, containing m — k letters A and k letters B. Although the conjecture is 
widely expected to be true, there are, by now, only partial results confirming it. Of course, 
it is true in the obvious cases of commuting A and B, for n = 1 and for k < 2. For n = 2, 
the statement follows since there is a common basis where A and B have nonnegative entries 
only [12j . Beyond that, positive results have been obtained for lower m; at present, the 
conjecture is proven for m < 13 [lOtlH]. This relied on two main ideas: First, generally, if the 
conjecture is given for some {ni,k), then it holds for any (m',k') with m' < m, k' < k and 
m' — k' < m — k [8\. Second, more specifically, Hagele [7] proposed to write Sm,k{A, B) — up 
to some cyclic permutations — as a sum of positive terms. Although not possible for (6, 3) 
and several other cases he was able to find such a decomposition for (7, 3), implying the 
BMV conjecture for m < 7. More refined methods |10J using computer algebra established 
the cases (14,4) and (14, 6), implying the conjecture for m < 13. Recently, it has been shown 
that the conjecture is always true for A: = 4 [3], implying it for A; = 3. Other results show 
that one may restrict oneself to the case of singular matrices A and B when proving the 
conjecture inductively [8]. Although the BMV conjecture is still open, it is known that the 
untraced coefficients Sm,kiA,B) need not be positive. The easiest example is 56,3(^,-6) for 
appropriate A and B; here, some single words may even have negative trace [9]. 

In the present paper we study a different side of the problem. Shifting the focus from 
(computer) algebra back to analysis, we are going to investigate the behaviour of the terms 
Sm,k{A, B) for large instead of small m. Our main result i^ 

Theorem 1.1 Let A and B be positive hermitian n x n matrices and A; G N. 
Then there is some mo S N, such that: 

AB = =^ tr Sm.ki^j B) = for any integer m ^ k ^ 0. 
AB / tic Sm,k{A, B) > for any integer m > mo. 

We are going to prove this theorem in two different ways — once using complex-analysis 
methods, once using combinatorics. In the latter case we also give a concrete estimate for 
rriQ. Let us summarize the main ideas of the proofs. Since the case AB = is trivial, we 
may assume AB ^ 0. Moreover, we may assume that A has unit norm0 

^Originally, in [1], a stronger conjecture has been stated: For any bounded- from-below self-adjoint opera- 
tors A and B and any eigenvector 95 of -B, the function {(^, e~'^^*^'(/3) is the Laplace transform of a positive 
measure whose support is contained in the convex hull of the spectrum of B. This conjecture, however, 
turned out to be wrong as seen by Froissart (see the notes added in proof in [1]; alternatively, see the example 
given in U). Then, BMV conjectured that, nevertheless, the statement remains valid for the trace. 

^Positivity of a hermitian matrix B means that (a;, Bx) > for all x £ C". In particular, is positive. 

■^If the dimensions of the matrices and 1 should be clear from the context, we may refrain from specifying 
them by writing On and 1„, respectively. 

''in the main text, we always consider the operator norm. Other norms will be discussed in the context of 
Euler-Lagrange equations in the appendix. 
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For the combinatorial proof observe that, if m increases, the k letters B are getting more 
and more sparsely distributed inside the words in Sm,k{^-, B). Indeed, most of the terms are 
of the form A^^BA'^ ■ ■ ■ BA^'=+^ with rather large i^. These words are approximated by the 
positive hermitian matrix {PaBPa)^ , where Pa is the hermitian projector limj^oo^*- The 
assertion follows unless tr {PaB)^ = tr {PaBPa)^ > vanishes. But, then A = © A' 
and B = Qn-l ® B' for some positive hermitian I x I matrices A' , B' with < Z < 77,, such 
that the proof follows inductively. 

For the complex-analytic proof, again by induction, we also may assume tr {PaB)^ > 0. 
Consider the series 

m=k ^ ' K=0 ^ ' 

Here, each c«; is a rational function vanishing at infinity and having poles only outside the 
closed unit disk. Now the proof follows, since the m-th Taylor coefficient of (1 — r)""'* is a 
polynomial in m of degree k — 1 with positive leading coefficient. 

Unfortunately, the dependence of mg on A and B is crucial for our proofs of the theorem. 
Therefore, the full BMV conjecture does not follow directly from the theorem above. Nev- 
ertheless, some (admittedly, simple) numerical simulations indicate further structures in the 
sequence of tr S'm^fe(^, S) for general k. To see them, we should first factor out the trivial 
dependencies. In fact, observe that otherwise this term (in general) diverges; we have, e.g., 
tr 5'^^fc(Kl, Al) = nK"^~^ X'' (Jj^y Thus, we will study the normalized quotient 

qmMA,B) :- tr5^,,(P||l,||i?||l)' 

as the BMV conjecture is now equivalent to qm,k{^^ B) >0 for all positive hermitian matrices 
A and B having norm 1. Since the theorem above tells us that qm,ki^j B) > for sufficiently 
large m, the BMV conjecture would now follow if one could establish 

Conjecture 1.2 Let A and B be positive hermitian n x n matrices with AB 7^ 0. 
Then, for any fixed A; G N, the sequence 

{Qm+k,k{^^B))men 

is decreasing. 

Despite to the mentioned numerical hints, we have not been able to prove this conjecture 
analytically. Nevertheless, we have been able to deduce further properties of qm,k{A,B) for 
large m and general k: 

Theorem 1.3 Let e > 0, let A and B be nonzero positive hermitian matrices, and let d be 
the dimension of the intersection of the eigenspaces of A and B w.r.t. their 
highest eigenvalues. Then there are mo, A;o > 0, such that 

qm,kiA, B) > — — £ for m > mo and k,m — k > 



n 



and 



qm,k{A, B) < — he for m > mo and k,m — k > ko. 



n 



In particular, tTSm,k{A,B) is strictly positive for all k, provided m is sufficiently large and 
the matrices A and B share a common eigenvector w.r.t. the respective maximal nonzero 
eigenvalue. 

Let us now sketch the idea of the proof for normalized A and B. If d > 0, we may 
decompose A and B into A' © 1^ and B' © 1^;, respectively, where A' and B' are positive 
hermitian with IIA'S'll < 1. Since 
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= tl Sm,k{A. , B ) -\- — tl Sjn,k{'^ni '^n) J 

we may assume d = 0, i.e., < 1- Moreover, by Theorem 11.11 we may assume that k 

and m — k are not too small. Now, the typical element among the Sm,k{A, B) terms contains 
a higher and higher number of subwords AB. The norm estimate ||A-B|| < 1 implies that 
1m,k{A,B)^ hence, the average contribution of a word to Sm,kiA,B) is getting arbitrarily 
small. 

Our paper is organized as follows: First, for completeness, we collect some simple properties 
of normalized positive hermitian matrices. Next, we study properties of certain power series 
whose coefficients are Hurwitz products or their traces. Then we use combinatorial methods 
to the calculate the number of words in A and B containing the subword AB a certain 
number of times, and derive estimates for these figures. In Section O we prove the theorems 
announced above. Finally, in Appendix El we derive the Euler-Lagrange equations in a 
slightly more abstract way than in [8j and extend these results to several norms. 

2 Some Algebra 

In this section we review the asymptotic behaviour of powers of positive hermitian matrices 
as well as of their products. Most importantly, we will recall that A^' for positive unit-norm 
matrices A always tends to the proj;ctoJl onto the highest eigenspace (i.e., for the eigen- 
value 1); powers of matrix products converge to projectors to common highest eigenspaces. 
Moreover, we derive some norm and trace estimates as well as some criteria for the product 
of two matrices to vanish. 

2.1 Power Limits 

Definition 2.1 For any nx n matrix A, define Pa '■= limj^oo^N if the limit exists. 
Obviously, we have Pi = 1 and Pq = 0. 

Lemma 2.1 Let A be any n x n matrix, such that Pa exists. Then we have: 

• Pa is idempotent; 
. APa = Pa = PaA; 

• Pax = X <^=^ Ax = x, where x G C""; 

• Pu-^AU exists for any invertible n x n matrix U, and it equals U~^PaU. 

The final statement implies that we often may restrict ourselves to the case of diagonal A, 
as long as we investigate Pa for hermitian A. 

Proof We have 

PaPa = (limi^oo (linij^oo ^^') = limi^oo{A' {ihuj^oo A^)) 
= limj_oo(limj^oo-4*+-'') = limj_^oo -P4 = Pa 

and, similarly, APa = Pa = PaA. Now, x = Pax implies 

Ax = A{Pax) = {APa)x = Pax = x. 
The remaining assertions are obvious. qed 

Definition 2.2 Let A be any n x n matrix. 

Then Ix{A) denotes its eigenspace in for the eigenvalue A. 

^Throughout the whole paper, any projector is assumed to be a hermitian projector. 
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Lemma 2.2 If A is hermitian with ||yl|| < 1 and if —1 is not in the spectrum of A, then 
Pa exists and is a projector. Moreover, im Pa = -^i(^)- 

Proof Consider A in diagonal form and use Lemma |2.1[ qed 

Lemma 2.3 Let Ai be an rii x rii matrix and A2 be an n2 x n2 matrix, such that Pai and 
Pa2 exist. Then Pai(bA2 exists and equals Pai © Pa2- 

Proof Use that (^1 © ^2)* = A\ © A^^. qed 
2.2 Phone Matrices 

Definition 2.3 A matrix A is called n-phone iff A is a positive, hermitian n x n matrix 
whose largest eigenvalue is 1. 

Recall that the norm of a positive hermitian matrix coincides with its largest eigenvalue. 
Lemma 2.4 Any nonzero projector is an n-phone matrix. 



Lemma 2.5 For any k € N+ and any n-phone matrix A, we have 

" ' '-A 



{A- Pa)'' = A'' -Pa and \\A'' - Pa\\ = \\A-Pa"'' 



Proof This follows inductively, using 

= ^A-Pa)H^-Pa) = {A'^-Pa){A-Pa) 
= A''+^ - PaA - A'' Pa + PaPa = A''^^ - Pa 

by Lemma r2.ll The norm equality now follows from HM'^H = ||M||'^ for any hermitian 
matrix M. qed 

Occasionally, we will decompose matrices into direct sums of matrices. When we simply state 
that some matrix B equals Bi(BB2, then we tacitly assume that there is some decomposition 
of C" into Xi®X2 = C^'"^-^! © C^'-^^^^ such that B\x, = Bi : Xi — > Xi. Furthermore, 
note that whenever we decompose several matrices into direct sums, we will always assume 
that all these matrices are decomposed w.r.t. one and the same decomposition of C". 

Lemma 2.6 Let A be an n-phone matrix and let Pa = '\-n-i © 0; for some < / < n. 

Then there is some < a < 1 and some l-phone matrix A' , such that A equals 
In-i © oiA' . Moreover, we have I > 0, unless A = 1„, and I < n. 

Proof • If / = n, then Pa = 0„, whence Ii{A) = by Lemma \T2\ i.e., II^H < 1. 

• If / = 0, we have Pa = In, i.e., hiA) = hiPA) = n and, therefore, A = In- 

• If < Z < n, then Pa = ( ) • A = ( ^ £f ) with positive hermitian matrices 
F (of size n — l) and H (of size I). From PaA = Pa, we derive F = \ and G = 0, 
whence A = In^i © H. By 

l„_i©Oi = Pa = Pl„_i(BH = Plr,^i®PH = ln~l®PH, 

we have Ph = 0;, whence \\H\\ < 1, again by Lemma l2.2[ Now, define a := 
and A' := a'^H (or A' = li H = 0^). qed 

Corollary 2.7 For any n-phone matrix A, we have \\A — Pa\\ < 1- 
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2.2.1 Shared Eigenspaces 

Lemma 2.8 Let Ai, . . . ,A]\f be n-phone matrices and let x £ C". Then 

\\Ai\f ■ ■ ■ Aix\\ = \\x\\ <^=^ AiX = X for all i = 1, . . . , N . 

Proof We may assume that x ^ 0. Moreover, the dhection is trivial. We now prove 
the =^ statement by induction. Let = 1 and denote shortly A := Ai. Then 
there is a unitary U, such that D := UAU* is diagonal. Setting y := Ux, we have 

||Dy|| = \\UAU*Ux\\ = \\Ax\\ = \\x\\ = \\Ux\\ = \\y\\. 

Writing D =: diag {di, . . . , dn) with Q < dj <1 and y =: {yi, . . . , yn)^ , we find that 
Dy = {diyi, dnynV , whence 

E-=i(i-^|)l%f = E-=il%f -E-=irf>.f = \\vr-\\Dy\? = 0. 

Since < < 1, we have (1 — = for all j. Consequently, 

dj = 1 or = for all j =^ {dj — = for all j 

=^ djyj = yj for all j 
=^ Dy = y. 

Now, Ax = U*DUU*y = U*Dy = U*y = x. 

Next, let > 1 and assume the assertion to be proven for — 1. We now have 

||x|| = P^^7V-i---^ia;|| < ||A^||p7V-i---^ix|| 

= II Aat-i • • • ylix|| < ||yliv-i|| • • • ll^ill ll^ll = ll^lli 
whence ||j47v-i • • • ^i2;|| = ||x||. By induction, AiX = x for all i < N. On the 
other hand, this implies II^atxH = || AatAat-i • • • Aix|| = ||x||. From the induction 
beginning, we get A]\[X = x as well. qed 

Corollary 2.9 For any n-phone matrices Ai, . . . , A^ we have 

h{Ai---AN) = /i(Ai)n...n/i(Ajv). 

Proof • Let x £ Ii{Ai ■ ■ ■ A^), i.e. Ai - ■ ■ A^x = x. Lemma ETH] implies AiX = x for all i. 

• Trivial. qed 

Corollary 2.10 For any n-phone matrices Ai, . . . , Aj\f we have 

hiAi---AN) / ^ \\Ai---An\\ = 1. 

Proof • If 11^1 • • • ^7v|| = 1) then there is some nonzero x £ C" with ||^i • • • ^Arx|| = ||x||. 

Lemma [2.81 implies that AiX = x for all i. This, of course, implies Ai - ■ ■ Anx = x. 

• If Ii{Ai ■ ■ ■ A]\r) / 0, then, by Corollarv 12.91 there is some nonzero x S C", such 
that AiX = X for all i. Now, Ai - ■ ■ A^x = x and ||^i • • • ^Ar|| = 1. qed 

Corollary 2.11 Let Ai, . . . jA^ be n-phone matrices. Then we have: 

-Pai - Ajv ^^^s*^ ^'^'^ ^'l^^ls 0- h{Ai • • • Aj\f) = 0. 

Proof • If Ii{Ai ■ ■ ■ Aj\[) = 0, then, by Corollary 12.101 ll^i " ' ' ^A^ll < 1) whence we have 
11(^1 • --ANyW < \\Ai ■ --AnW^ for j oo. Consequently, Pa-^-An = 0- 

• If Ii{Ai ■ ■ ■ A]\f) / 0, then, by Corollary 12.91 there is some nonzero x S C", such 
that AiX = X for all i. This means, Ai - ■ ■ A^x = x and thus Pai---Ai^x = x. 

qed 
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Lemma 2.12 For any n-phone matrices Ai, . . . , A]\[ we have: 

1. There are n' x n' matrices A'^, . . . , A'j^, such that for i = 1, . . . , N 

• Ai = A[® h, 

• each j4. is positive hermitian; 

• \\A[---A'^\\<1. 

Here, / := dim/i(j4i • • - An) and n' := n — I. 

2. Paj^-An exists and is the projector to Ii{Ai ■ ■ ■ An). 

Proof Denote Ii{Ai ■ ■ ■ An) C C" shortly by X. By Corollary EJl each Ai is the identit 
when restricted to X. Since each A^ is hermitian, X-^ is preserved by each A^ 
Hence, we may decompose each Ai into Ix © A^ according to C" = X © X^. Here, 
A'- is a positive, hermitian operator on X-^. (W.l.o.g., we may assume that A'- is an 
n' X n' matrix with n' := n — dimX.) If \\A[ ■ ■ ■ A'j^\\ was 1, then 

1 ~ 11^1 ■ ■ ■ ^Afll — 11^1 II ■ ■ ■ II^Afll — ^1 

and each A'^ would be n-phone. Since, however, by construction and by Corollary 
EHl = Ii{A[) n . . . n hiA'j^) = Ii{A[ ■ ■ ■ A'j^), we have Pa\---A'j^ = 0> as shown in 
Corollarv I2.11[ Consequently, by Corollarv 12.101 M'l " " " ^ivll / ^- Obviously, we 
have Pi^ = Ix, such that, by Lemma Pai---An = P{ix®A[)--iix(BA'^) exists and 
equals Pi^ © Pa[ - a'^ = Ix © Ox-l- It is, of course, hermitian. qed 

Corollary 2.13 Ii{Ai ■ ■ ■ An) = Ii{Pai-An) any n-phone matrices Ai - ■ ■ An- 
2.2.2 Norms and Traces 

Lemma 2.14 Let A,B be n-phone matrices. Then < for all i > j. 

Proof Let D be the n-phone matrix with B = D^. Then 

\\AB'A\\ = \\{D'A)*{D'A)\\ = \\D'Af < \\D'~^f\\D^Af = \\AB^A\\. 

qed 

Corollary 2.15 Let A,B be n-phone matrices. 

Then ABA = implies AB'^A = for any z G N+. 

Proof We have = \\ABA\\ > p5M|| > by Lemma EH qed 

Lemma 2.16 For any n-phone matrices A and B, we have 

\\BABf+^ < ||(^P^)'=|| < \\BABf-\ 
If B is even a projector P, then 

\\PAPf < ||(^P)''|| < \\PAPf-\ 

Proof Since BAB = B*AB is hermitian and positive, we have ||(i?^-B)'^|| = ||B745||'^ for 
any k G'N. Now observe that 

\\BABf+^ = \\{BAB)''+'^\\ = \\B{AB^)''AB\\ 

< \\{AB^)''\\ = \\AB{BAB)''-'^B\\ 

< \\{BAB)''-^\\ = \\BAB\\''-\ 

since ||j4|| = 1 = ||i?|| and, similarly. 



^Let e and x G X. Then {x,AtX-'-) = {A*x,x^) = {AiX,x^) = {x,x-'-) = 0, hence A^x^ G X^ 



7 



llP^Pf = ||(P^P)^|| = WPiAPfW < WiAPfW = \\AP{AP)''-^\\ 
< \\P{AP)''-^\\ = \\iPAP)''~^\\ = \\PAPf-\ 

since P^ = P and ||P|| = 1. qed 

Proposition 2.17 Let A and B be n-phone matrices, and let k E N+. Then 

AB = tvAB = tr (AB)'' = ABA = 0. 

Proof Let C be an n-phone matrix with A = C'^. 

• First of ah, let tr (AB)^ = for some k E N-j.. Since 

ir (AB)'' = ti C{CBCf-^CB = ti (CBC)'' 

and since CBC is positive hermitian, we hav^ CBC = and ABA = 0. 

• Next, ABA = imphes {BA)*BA = AB'^A = by Corollary [2T5l whence 
\\BAf = \\{BA)*BA\\ = 0, implying BA = Q and AB = 0. 

• Finally, of course, AB = implies tr {AB)'^ = 0. qed 



2.2.3 Splitting 

Lemma 2.18 Let A and B be n-phone matrices. Then AB = iff there is some < / < n, 
some /-phone matrix A' and some (n — /)-phone matrix B' , such that 

A = A' ®Qi and B = 0„_i B' . 



Note again, the splitting above means that there is a basis of C", such that A and B can be 
simultaneously splitted in the way given above. 

Proof If A and B can be split in the given way, then AB obviously vanishes. The other 
way round, AB = implies BA = 0, hence AB = BA, whence A and B can be 
diagonalized simultaneously. Now, the statement is trivial. qed 

Lemma 2.19 Let k E N+, and let A and B be n-phone matrices. 

Then we have tr (PaB)'' = iff there are 0<a<l, 0<l<n and /-phone 
matrices A' and B' with 

A = In-i © aA' and B = Q^-i ® B' . 



Proof By Proposition 12.171 ii {PaBY = is equivalent to PaB = 0. Analogously to the 
proof of Lemma l2.18t we see that, for PaB = 0, there is a decomposition 

Pa = In-i © 0/ and B = 0„_z © B' 

for some /-phone matrix B'. Since Pa 7^ 1 by PaB = 0, we have < / < n. Now 
the implication follows from Lemma |2.6[ The other direction is trivial. qed 



2.3 Hurwitz Products 

In this subsection, for completeness, we list several properties of Hurwitz products and their 
traces. The proofs are simple and therefore omitted. They may also be found in 

Definition 2.4 Let m and k be integers, and let A and B he n x n matrices. 

Then the Hurwitz product Sm,k{A, B) is the sum of all matrix products 
containing exactly m — k factors A and k factors B. 

^tr = implies _D = for positive hermitian matrices D. 
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For definiteness, we assume the Hurwitz product to be zero if m, k, or m — k is negative. 

Lemma 2.20 The Hurwitz product of any two hermitian n x n matrices is hermitian. 
Consequently, its trace is always real. 

Lemma 2.21 For any m, € N and any hermitian n x n matrices A and B, we have 

Sm,k{A.,B) = ASm-i,k{A, B) + BSm-l,k-l{A, B), 
Sm,k{A,B) = Sm-l,k{A, B)A + Sm~l,k-l{A, B)B. 

Lemma 2.22 For any m, S N and any hermitian n x n matrices A and B, we have 

{m - k) tr SmM"^^ ^) = rntrASm-i,kiA,B), 

ktrSm,kiA,B) = mtr BSm-i,k-iiA,B). 

Lemma 2.23 Let Ai and Bi be hermitian rii x rii matrices with rij G N for i = 1,2. Then 

Sm,k{Ai(BaA2,Bi®(3B2) = 5i) + a'"-'=/3'= S^,fc(^2, ^2) 

for all m,k G N+ with m > k and a, /3 G C. 

3 Some Complex Analysis 

To prove Theorem 11.11 using complex-analytic methods we will need to study the behaviour 
of 

tr [5(1-tA)-1]^ 

for n-phone matrices A and B. In this section, we provide the necessary statements. 

Lemma 3.1 Let A and B be n-phone matrices and let A; G N. Then we have 

00 

{1-tA)-'[B{1-tA)-']' = Y,T^S^+kAA,B) 

m=0 

for all r G C with |t| < 1. 
Proof Since ||^|| = 1 and |t| < 1, 

m=0 m=0 

converges absolutely and gives the assertion for /c = 0. Inductively, we have 

00 00 
(1 - tA)-^ [B{1 - tA)-^] = ^ T™ Sm+k,k{A, B)B 

m=0 m'=0 

00 

= {Sm+k+i,k+i{A, B) - Sm+k,k+i{A, B)A) A^ 

m,m,'=0 

00 00 

= Y r"'^"''Sm+k+i,k+i{A,B)A^' - Y r^^'^'Sm+k+i,k+iiA,B)A^' 

m,m'=0 m=0, m'=l 

00 

= T"^Sm+k+l,k+l{A, B), 

m=0 

where we used 

Sm+l,k+liA, B) = Sm,k+iiA,B)A + Sm,kiA,B)B 
from Lemma [2.21l in the second step and Sk^k+iiA, i?) = in the third one. qed 



Corollary 3.2 Let A and B be n-phone matrices and let k £ N+. Then we have 

k 

m + k 



oo 



i.[B{l-TAr^]' = ^r"^-^tTS^^,MA,B) 



m=0 

for all r G C with Irl < 1. 



Proof Use the relation k tr Sm,kiA, B) = m tr i?5m-i,fc-i(j4, B) in Lemma 12.221 to derive 
ix[B{l-TA)-^]^ = tiB{l-TA)-^[B{l-TA)-^]^'^ 



oo 



y^T^'tiBS^+k-iM-M.B) = Vr™— ^tr5„+fc,fc(Ai?). 



m + k 

m=0 771=0 



qed 



Lemma 3.3 Let k he a. positive integer, and let be rational holomorphic functions for 
K = 0, . . . ,k — 1. Assume that there is a real r > 1, such that none of the Cf^ 
has a pole for |r| < r. Finally, let G C have positive real part and define 
the analytic function / and its expansion coefficients fm by 

m=0 ^ ^ K=0 ^ ' 

Then there is an mo G N, such that Re fm> ^ for all m > niQ. 
Proof Write / as a Laurent series 

around 1 for appropriate G C and holomorphic F. Then F is again rational and 
has no pole for |r| < r. Consequently, the norm of the m-th Taylor coefficient of 
-F(r) can be estimated by ^ for some constant C > 0. As, by Lemma IC.21 the 
m-th Taylor coefficient of 



1 X ^ / m + K — 1 



(1 — t)'^ ^ \ m 

^ ' m=0 ^ 



is a polynomial in m of degree k — 1 with leading coefficient ■, the assertion is 

obvious since Re > 0. qed 



4 Some Combinatorics 

The ultimate goal of this article is to derive asymptotic properties of in Srn,k{,A^B). Recall 
that 5m,fc(^, -B) equals the sum of all products of matrices where m — k factors equal A and 
k factors equal B. The trace of such a single product significantly depends on its "factor 
pattern". For instance, if the substring AB appears I times in the matrix product, then the 
trace of the full product cannot exceed n||^i?||'. To finally estimate the sum of all these 
product traces, we need estimates how frequently this pattern appears. This now is a purely 
combinatorial problem for words in two letters. To avoid confusion we will denote the letters 
by a and 6, and return to A and B only later. Let, moreover, < < m be integers and 
denote the set of all words containing exactly m — k letters a and k letters h by Wm,k- 
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4.1 Counting 



Proposition 4.1 Denote by Cm,k,s ^ y^m,k the set of words containing exactly s times tlie 
subword ab. Then we liave 



m,ky 



m — k\ fk 



and 



in,k 



Here, we used tlie convention tliat Q) = if j > i. 

Proof Let w G Cm,k,s be a word with exactly s subwords ab. Then 

for appropriate z^, > 1, i = 1, . . . , s, and jo, is+i > with ii + . . . + i^+i = m — k 
and jo + . . . + js = /c. Obviously, it is sufficient to prove that there are exactly (^) 
ways to write /c as a sum jo + ji + • • • + js of s + 1 integers with jo > and jt > 1. 
In fact, there are (^) possibilities to choose s elements Ji < . . . < Jg out of the 
k numbers 1, . . . ,k. Letting jo := Ji — 1 and j^ := J^+i — for < t < s and 
js ■= k + 1 — Js gives such a decomposition jo + ji + • • • + js of k. Since, the other 
way round, each such decomposition can be obtained by such J^, we get the proof. 
The second assertion is clear. qed 



Lemma 4.2 1. No word in Cm,k,k contains the subword 66, i.e., any word in Cm,k,k can be 
written as a^^ba^'^ ■ ■ ■ a^''ba^^+'^ with «t > and ik+i > 0. 
2. Denote by 'Dm,k,L ^ Cm,k,k the set of those words a^^ba}^ ■ ■ -a'^^fea^^+i as 
above with L and ik+i > L for some integer L > 0. Then 

\T^m,k,L\ = \Cm-(k+l)L,k,k\- 



Proof 1. This follows directly from the proof of Proposition 14. 1[ In fact, let 

w = V'^a'^V^ ■■■a'^b>'^a"^+^ G Cm,k,k- 

Since jo + • • • + jfc = A; and ji, . . . , j^ > 0, we have jo = and ji = . . . = j^ = 1. 
2. One easily checks that 

^ • ^m-{k+l)L,k,k ^ '^m,k,L ^ Cm,k,k- 

a}^ba^^ ■ ■ ■ a^^ba^k+^ i — > a'i+^6a'2+^ . . . a'*+^5a*'=+i+-^ 
is a bijection. qed 



4.2 Estimates 

We will need estimates on how the number of words changes in the event of having a fixed 
amount of letters less and how often there are subwords ab. In the first simple lemma, we 
will see that the (relative) decrease of the word number while dropping a finite number of 
letters a is arbitrarily small provided we had started with a occurring sufficiently often. In 
the second lemma, we show that — again for a occurring sufficiently often, i.e., for large m — 
the (relative) number of words containing less than k subwords ab can be made arbitrarily 
small. Or, in other words, if one of the k letters b appears then it appears "lonely", i.e., b^ 
or higher powers typically do not appear. 



11 



Lemma 4.3 For < e < 1, positive integers L and m with m > L{1 + -), we have 

m — L\ , fm 



Proof Use 



r\ / \ L-l , . \ L-1 / 

m — L\ I m\ -TT / K \ I m\ -i-r £\ Im 



k J \ k J \ m — j J \ k J -^^ \ LJ \ k 



qed 



Lemma 4.4 Let < e < 1 < 5 and 
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m > h 25* — 1 and k,m — k > S. 

£ 



Then 

5-1 



E 

s=0 



m — k\ fk\ f m — k\ fk 



s J \ s J ^ V S J \S 



Proof Observe that for 0<s<S<k,m — k and for m as in the lemma 



s 



2 
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— < — < — + 2S'-l-(2s-l) < m-2s + l 

e e £ 

< {{m- k) - s){k - s) +m-2s + 1 = {{m - k) - s + l){k - s + 1). 
Using the abbreviation dg := {"^J^) (^) for all s G N, one immediately checks that 

^''^ ~ {m-k-s + l){k-s + l) 
As just seen above, the prefactor is always smaller than e < 1, whence we get 
S-l S-l ^3 

s=0 s=0 

qed 

5 Proofs of the Main Theorems 

5.1 Reduction of Theorem 11.11 to the Case PaB 7^ 

Proposition 5.1 If the assertions of Theorem 1 1 . 1 1 hold for m > k > and for any n-phone 
matrices A and B with PaB 7^ 0, then Theorem II .11 is valid in toto. 

Proof For the proof of Theorem II. H we may, of com'se, assume that A and B are n-phone 
matrices. Moreover, we may assume that m > k > 0, that n > 1 and that AB 7^ 0, 
as the other cases are trivial. 

Assume now that PaB = 0. Then, by Lemma 12.191 we find some < a < 1 and 
some /-phone matrices A' and B' with < Z < n, such that 

A = In-i © a A' and B = 0„_/ © B' . 

Since ^ AB = 0„_/ © aA'B', we have q 7^ and A' B' ^ 0;. Together with 

Sm,k{A,B) = Sm,ki^n-l,On-l) + Oi"^"'' Sm,kiA' ,B') = a"^^^ Sm,kiA' , B') 

by Lemma 12.231 and m > k > 0, this implies the assertion by induction. qed 
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5.2 Complex-Analytic Proof of Theorem 11.11 

First we prove our Main Theorem by means of complex analysis without focussing on concrete 
estimates. 

Proof Theorem 11.11 

Let A and B be n-phone matrices. Using Lemmata 12. II and 12.51 we have 



m=0 m=0 

and 

oo 

{1-Pa){1-tA)-' = Y^il - Pa) t"' A"' 

oo oo 
m=0 m=0 

Observe that (1 — Pa){^ — tA)~^ is rational and (up to the removable discontinuity 
at 1) analytic for |r| < — PaH""*^, whereas — PaH""*^ is strictly larger than 1. 
From 

B{1-tA)-' = B{PA + {l-PA)){l-TAy' 

= -^BPa + B{1-Pa){1-tA)-\ 
i — r 

we derive 

whereas each is a rational function vanishing at infinity and having poles only 
for 

, , 1 

\\A-Pa\\ 

Now the assertion follows from Lemma 13.31 together with Corollarv l3.2t In fact, we 
may assume PaB ^ 0, hence trP^i? > 0, as well as /c > by Proposition 15.11 and 
we know that any Hurwitz product trace is real by Lemma |2.2(JI qed 

5.3 Combinatorial Proof of Theorem 11.11 

There are two main steps in the study of the asymptotics of tv Sm,kiA, B) for growing m 
while k is fixed: First, we estimate how fast the products A'-^B • • • yl'^'S do approach {PaB)^ 
depending on a lower bound L to all /j. Second, the longer the words are (i.e., for growing 
m), all other words (i.e., those with < L or having substrings 6^) get less frequent for fixed 
L. This allows us to estimate how fast tr 5m,fc(^7 -B)/tr 5m.A;(l, 1) approaches ii {PaBY /n 
and, finally, to prove Theorem ll.il 



Lemma 5.2 For any n-phone matrices A and B and for any integers h, . . . ,lk > L > 0, 
we have 

k k 

l^^YlA^'B ~Y[PaB < k\\A- PaW^ WA^Bf-^. 

1=1 i=l 

Proof Observe that for any n x n matrices Xi, . . . , Xk and X, we have (see Lemma lC.ip 
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i=l 

Now, Lemma 12.51 implies 

k k k 

\\^Y[A^^B -WPaB < ^\\{PaBY-^\\\\{A^^ - Pa)B\\\\A^^+^B\\---\\A^''B\ 

i=l i=l i=l 

k 



< llP^^ir"^ \\A - PaW^' \\A^'+^B\\ ■ ■ ■ WA^^Bl 
1=1 

k 



i=l 

= k\\A- PaW^ \\A^B\\''-\ 

qed 

In Section HI we studied words in the two letters a and b. We now define W to be the 
homomorphism from Wm,k to the n x n matrices, whereas W{a) := A and W{b) := B. It is 
now clear that, e.g., Sm,kiA,B) = Z]^„gw„,fc W{w). 

Proposition 5.3 Let A and B be n-phone matrices, and let A; € N and e € (0, 1) be fixed. 
Choose now some L S N+, such that 

k\\A- Pa\\^ \\A^Bf-^ < e. 

Then, for any m E N with 

k 



m 



> + (k + kL + L), 

we have 

tTSm,k{A,B) lT{PABf ^ / triPAB)'' ^ ^ 



tr5m,fc(l,l) n 



n 



Observe that tr (PaB) is always nonnegative. 

Proof First observe, that \\A — Pa\\ < 1 by Corollarv 12.71 whence there exists such an L. 
Next, observe that 

\y^m,k\T^m,k,L\ = |VVm,fc| " |^m,,fc,L| 



- \Cm-{k+l)L,k,k\ (by Lemma 

m\ fm- {k + l)L -k\ fk\ , , , 

k)~[ k )[k) (by Proposition iH) 

/ 7Tl\ ( Tn\ 

< ( ^ j - (1 - e) L j (by LemmaESl) 

= e|>Vm,fc|- 

since m > ((A; + 1)L + A;) (l + ^) by assumption. Third, observe that for w € ^?m,fc,L) 
we have 

\tr W (w) -ti {PaB)''\ < nk\\A- PaW^ \\A^Bf-^ < ne 
by Lemma 15.21 and |trC| < ^^||C|| for any matrix C, whence 
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< 



< 



tr 5^,^(1,1) n 

s-^ tTW{w) -tT{PAB)'' 

trW{w) -ti{PAB f 



E 



ne + 



< 

using 



\'^m,k,L\ 

n\Wm,k\ 

tr (PaB)'' 

n 



tr ^^^^(l, 1) 

|VVrn,fc| — |^m,fc,L 



m,k I 



(n + |tr(P^B)*^|) 



trT^(u') - tr {PaBY 



tr^m ,fc(l, 1) 



E 1 



- tr S"™ 1) 

n 



qed 



Theorem 11.11 is now a corollary: 



Proof Theorem 11.11 

By Proposition 15.11 we may assume that fc > and that A and B are n-phone 
matrices with PaB 7^ or, equivalently, with tr(P^i?)'^ 7^ by Proposition 12.171 
Then there are L > and e £ (0, 1) with 

tr {PaB)'' 



k\\A- PaW^ WA^Bf-'^ < e < 



Now, since < tr (PaB)'' < n, we have 

tlSmA^^B) ^ tl{PABY (iT{PABf 



tr 5^,^(1, 1) 



> 



> 



n 



3n 



+ 2 e 



n 



ti {PaBY ftv{PAB 



n 



+ 2 



n 



1 tr (PaB) 



n 



tr (PaB)' 



n 



trjPABY 
3n 

> 0, 



provided 



m > mo 



(1 + {k + kL + L) 



qed 



Remark The proof of Theorem 11.11 above provides us with an explicit estimate for the 
value of rriQ. If A is not a projector and PaB 7^ 0, we have tr 5m,fc(^> i?) > for 
all m > tuq with 



mo 



(1 + A;) 1 + 



3kn 



2 + 



lntr(P4S)'' -ln3A;n\ 



InlM -PaI 



tr (PaB)''. 

If A is a projector and AB ^ 0, then tr S„i,k{^, B) > for all m > mo with 

3kn 



mo 



1 + 



2k) 



1 + 



tr (AB)''. 

For PaB = 0, use the decompositions of A and B as in the proof of Proposition 
15.11 and then use the expressions above with A and B replaced by A' and B', 
respectively. (If again Pa'B' = 0, proceed iteratively.) Of course, the estimates 
above need not be optimal; if the BMV conjecture was true, mo would probably 
be k unless AB = 0. 
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5.4 Growing m and Not-too-small k 



If ab appears S times in a word in Wm,k, then the corresponding matrix product has at 
most norm P^ll^. For growing m, the typical number of alternations between o and b in 
a word indeed increases; in particular, it passes the threshold S sooner or later. Therefore, 
the normalized trace of Sm,k{A, B) can be estimated by up to some e. 

Proposition 5.4 Let < e < 1 < S" for some 5" € N and 
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m > h 25 — 1 and k,m — k > S. 

e 



Then 



tr5^,fc(A5) 



tr5m,fc(l, 1) 



< e + \\AB\ 



Proof First observe that |trl^(it;)| < n||VF(z«)|| < for any w G 'Dm,k,s with s > S. 

Now, we simply decompose all elements of Wm,k into two sets: one consisting of 
all elements containing less then S subwords ab and the other one consisting of all 
elements with at least S subwords ab. We get 

\tTSm,kiA,B)\ < iti'^Hi + E E itrw^Hi 

< \Cm,k,s\ 

S<S S>S W£Cm,k.s 



n+ > ^ n II^-BII-^ 

£ ,k.s 

-'m,k,S\ n + \ Wm,k\ n\\AB\\'' (by Lemma O]) 



< n|W^,fc|(e+ \\ABf) 

< trS„,fc(l,l) (e+ ). 



qed 



5.5 Asymptotics for Growing m and General 

Since Ii{A) f] Ii{B) = Ii{AB), Theorem 11.31 follows immediately from 

Theorem 5.5 For any n-phone matrices A and B, and for any < e < 1, there are some 
niQ G N and some /cq G N, such that for all m > mo 

tTSm,k{A,B) dim h{AB) , n n ^ ^ 

> e for all < K < m 



n 



and 



— — — < h £ tor all /co < K < m — Ko • 

tr6m,fc(l,l) n 



Proof First of all, let us find /cq and ttiq, such that 

ti Sm,k{A^B) dim h{AB) 



< e 



trS'm,fc(l,l) n 
for all m > ttiq and k^ < k < m — k^. 

• Assume first ||^-B|| < 1, i.e., Ii{AB) = by Corollarv l2.10i Choose some positive 
integer k^, such that 
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< 



and some mg G N, such that 

> + 2ko - 1. 

Now, Proposition 15.41 imphes that 

trS'm,fc(l, 1) 
for ah m > uiq and all ko < k < m — ko. 

Assume now = 1. According to Lemma [2.12l we may decompose A and B 

into A = A' ell and B = B' ® li with \\A'B'\\ < 1 for / := dimIi{AB). Using 
Lemma |2.23|, we have 

tT:Sm,k{A,B) = tT:Sm,k{A',B') +tlSm,ki'i-lAl), 

and, therefore, 

tvSm,k{A,B) I _ I ti S^,k{A',B') 



trS'm,fc(l„, 1„) n n tr ^^.^(Ij, 1;) ' 
If II A' II II II = 1, then A' and B' are /-phone matrices with ||A'i3'|| < 1, for 
that the result has been estabhshed above. 

If P'llll^'ll < 1, then choose ko G N, such that (|| A'|| H^'H)''" < e. Then, by 
Lemma l2.23| we have 



tiSm,k{A,B) I 



tr 5"^ In) ^ 

for any m,k gN with m — k > ko and k > ko 
Now, let us finish the proof by showing that 

tr Sm,kiA,B) dim h{AB) 



< 1 M'lr-'^llB'f < e 
n 



> 



trS'„_fc(l,l) n 
for all m > mo with an appropriate mo and for k < ko or k > m — ko. 

• Assume again first that ||^-B|| < 1. Now, according to Theorem 11.11 for each 
k € N, there is some integer m'Q{k), such that 

trS^,kiA,B) > and tr Sm,m-k{A, B) = tTS^,k{B,A) > 

for all m > m'o{k). Now, simply define 

nio := max {mo(A;), ttiq}, 

A;<fco 

and we have the desired assertion. 

• If ||^-B|| = 1, we decompose A and B as above into A = A' ®\i and B = B' (Bli 
with ||yl'i?'|| < 1. Again using Lemma 12. 23^ the assertion follows as in the 
previous case. qed 
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Appendix 



A Euler-Lagrange Equations 

In the main body of the article, we have always used the operator norm for matrices and 
reduced our investigations typically to normalized matrices. In fact, this has been justified 
by the homogeneity of Sm,ki'j ■)• Nevertheless, there is a full range of other possible norms 
that can be taken to normalize the matrices. In [8], e.g., the Frobenius norm has been used 
to derive the Euler-Lagrange equations of the BMV conjecture. They yielded, among others, 
relations between and ASm-i,k{A, B) in any point where tr Sm,k is minimal or maximal. 
In this appendix we are going to rederive these relations in a slightly more abstract way and 
extend them to other norms. 

For that purpose, we choose the following Schatten p-norm^ 

\\A\\p := \/trAP and \\A\\oo := \\A\\ 

for p > 1 and for positive hermitian A. One immediately sees that ||^||oo = linip^oo ll^llp- 
Let us now fix some p S [l,oo]. Moreover, to avoid cumbersome notation, we let n-phone 
matrices be positive hermitian matrices having p-norm 1 (instead of to be of operator norm 
1 as in the main text). Next, observe that for any matrix- valued functions / and g on M, we 
have 



tr{f + tgr = mtr ' (/ + tg) 

ax V ax 

and, by comparison of coefficients, 

tr5'^,fc(/,ff) = rnti {f Sm-i,k{f^s) + g' Sm-i,k-i{f,g))- 

Here, we abbreviate /' := etc. We now use two different types of functions for / and g: 
on the one hand, we keep the eigenvalues by conjugation with unitary matrices, on the other 
hand, we modify them by multiplication with appropriate commuting matrices. Namely, let 
first 

f(x) := e-^^^e^^ for C £ su(n), 

i.e., C* = -C and trC = 0. Then /(O) = A and /'(O) = \A,C\. Moreover, /(x) is n-phone 
for any x and any n-phone A. If now (j4, is an extremal point for tr 5m,fe among the 
positive matrices with unit p-norm, then we have for all C S 5u(n) 

= tr5;„,,(/,i?)U.=o = mtr([A,C]5^_i,fc(A5)). 

Since A and 5^-1,^(^,-6) are hermitian [S], we get [5m,-i,fc(^, -B), ^] = from Lemma fB.ll 
Now, secondly, we consider 

/(a;) := forCGs((n). 

Note that / may fail to be differentiable at x = for p = oo. In fact, let Eij be the matrix 
having entry 1 at position (i, j) and zeros elsewhere. Consider A := E\\ + E22 and C := En. 
Then ||^e^'^|| equals e^ for x > and 1 for x < 0, which is obviously not differentiable. 



'Note, that the Schatten p-norm is actually defined [2] by 



where Sj{A), j = 1, . . . ,n, are the singular values of A. For positive hermitian matrices, however, our notion 
coincides with that definition. As we are interested in the case of positivity only, we may sloppily re-use the 
notion p-"norm" for our case. In fact, our definition does not give a norm on the linear space of all n x n 
matrices. 
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In general, the problem arises if the maximal eigenvalue of A is of multiplicity 2 or higher. 
Therefore, for the moment, we assume p to be finite. One easiljifl checks that 

/'(O) = — ^-^{AC ti AP - A trAPC). 

Of course, a priori, it is not clear that f{x) is positive and hermitian, even for small x. But, 
if U is some unitary matrix, such that U*AU (and C/*5m-i,fc(^, -B)C/) is diagonal, then / 
is positive and hermitian for any C = UDU* with D being diagonal and real. In fact, the 
product of diagonal positive and hermitian matrices has these properties again. If now A 
and B are nonzero and again extremal for tr Sm,k among n-phone matrices, then 

= ^^tr5;,fc(/,i?)U=o = tr {{AC tiAP- A tiAPC)S^.i,k{A,B)) 
m ' ] 

= tr {Sm-i,k {A, B)A tr AP - AP tr Sm-iA^, B)A) C . 

Since S'm-i,fe(^, B) and A commute as seen above and are hermitian, we get 

Sm-i,k{A,B)AtiAP = APtj:S^.i,k{A,B)A 

from Lemma lB.21 Similarly, we can derive 

Srn^i,k-i{A,B)BtrBP = BP ii Sm-i,k~M,B)B. 

Altogether we have 

Proposition A.l If < < m and if tr Sm,k is extremal at (A, B) for the positive hermi- 
tian matrices having unit p-norm with 1 <p< oo, then 

S^^iAA,B)A = APtTS^^i,k{A,B)A 
S^^i,k~i{A,B)B = BPtvSrr,^^,k-i{A,B)B . 

The case p = 2 has already been derived by Miliar in 0. There, the norm equals the Fro- 
benius norm. The case p = oo, i.e., the supnorm case, can be dealt with as for p < oo as 
far as we derive that A and Sm~i,k{A, B) commute. Assuming now, for simplicity, that A 
and Sm-i,k{A, B) are diagonal and \\A\\ = 1, we see that f{x) := Ae^'^ is (at least for small 
|x|) n-phone — provided C is diagonal with Ca = for An = 1. Then /'(O) = AC implying 
tv AC Sm-i,k{A, B) = 0, whence the {i,i) components of Sm-i,k{A, B)A vanish if An ^ 1. If 
Pa has a single nonzero entry, we immediately get S'm-i,fe(^, B^A = P^tr Sm-i,k{A, B)A. In 
the other case, however, we run into the non-differentiability problem as above. At present, 
we are not able to solve this problem. 
Nevertheless, we have 

Corollary A. 2 If < A; < m and if tr Sm^k is extremal at {A, B) for the positive hermitian 
matrices having unit p-norm with 1 < p < oo, then 

S^AA,B) = {m-k)AP + kBP 

m 

The same is true for p = oo, provided P4 and Pg have rank 1. 
Proof Use the properties of Hurwitz products listed in Subsection 12.31 qed 



Observe that 

Pe^^||;(0) = i|lA|l^''(|lAe^^|l^)'(0) = ip|l^-(tr(Ae-^r)'(0) = UWl^ tr A^C 



and, therefore. 



1 i-P.. 1 



/'(O) = jij^{AC\\A\\,-A\\A\\l-^tvA^C) = —±-^(AC\\A\\;-AtrA^C). 
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B Lie Algebra Relations 



Lemma B.l If A and S are hermitian matrices, fulfilling tr [^4, C]S = for all C G su(n), 
then A and 5 commute. 

Proof Since A and S are hermitian, we have [A, S]* = —[A, S] and, anyway, tr [A, S] = 0. 

Therefore, [A, S] G su(n). Moreover, tr [S, A]C = tr [A, C]S vanishes by assumption 
for any C G su(n). Since su(n) is semisimple, the Killing form {X, Y) := ^ tr {XY) 
on su(n) is non-degenerate, giving [S, A] = 0. qed 

Lemma B.2 If A and S are hermitian matrices that are diagonal after conjugation with 
U and fulfill tr {{SA tr A^ - A^ tr SA)UDU*) = for all diagonal matrices 

D, then SA tiA^ = A^ tiSA. 

Proof If A and S arc already diagonal, then the assertion is trivial. In fact, letting D be 
the matrix having just a single nonzero entry at position the trace equation 

above means that the (i, i) component of {SA tr A^ — A^ tr SA) vanishes. Since the 
off-diagonal elements are zero anyway, we get the assertion. 
In the general case observe that 

tr {U*SU U*AU tv{U*AU)^ - {U* AU)^ tr U* SU U*AU)D 

= ii {S A ii A^ - A^ iiSA)UDU* 

reduces this case to the first one. qed 

C Simple, But Useful Identities 

Lemma C.l For any n x n matrices Xi and X, we have 

k 

X\ - • • Xk = X^ + ^ [Xi — X) Xi^i ■ ■ ■ Xk. 

i=l 

Proof For A; = 1, we have Xi = X^ + X^ {Xi — X). For A; > 1, we have by induction 

k 

Xi - ■ ■ Xfe+i = X'^Xk^i + ^2 {Xi — X) Xi^i ■ ■ ■ 

i=l 



= X''{X + {Xk+i -X))+Y, X'-^ {Xi - X) Xi+^ ■ ■ ■ XkXk+1 

i=l 

k+1 

= X'^X + ^ X* ^ {Xi — X) Xi^i ■ ■ ■ XkXkj^i. 



i=l 



qed 



Lemma C.2 For any natural number k and for all |r| < 1, we have 



1 fk + m 



(l_7-)fc+i 



Proof The statement is clear for A; = 0. Now, we have inductively 

k _ / 1 y _ ^ fk-l + m\ 

(l_^)fe+i - 1(1 -r)W ~ rn J 

\ / \ / m=l ^ ^ 

m=o ^ ^ m=o ^ ^ qed 
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